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An  ultraviolet  beam  (337.1  nm)  from  a  pulsed  nitrogen 
gas  laser  fired  through  a  fluid  containing  a  photochromic 
dye  produces  opaque  traces  in  the  fluid  along  the  path  of 
the  laser  beam.  Subsequent  movement  of  the  fluid  deforms 
the  opaque  traces  and  produces  a  record  of  the  fluid  flow. 

A  high  powered  laser  has  been  designed,  built,  and  tested. 
Velocity  distributions  have  been  obtained  in  slowly  forming 
droplets  of  solvent  in  an  immiscible  ambient  fluid. 
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I.  INTRODUCTION 


Valuable  insight  into  the  hydrodynamics  of  fluid  flow 
can  be  obtained  by  "seeding”  the  flow  with  some  tracer 
material  and  by  analyzing  the  subsequent  movement  of  the 
tracer.  Aluminum  flakes  have  been  used  to  demonstrate  sur¬ 
face  flows  in  liquids  [1].  Food  dyes  are  often  used  in 
water  and  smoke  streaks  can  provide  histories  of  air  flow. 
These  flow  visualization  techniques  provide  excellent  quanti¬ 
tative  measurements  of  fluid  flow  phenomenon,  and,  through 
the  use  of  photography,  qualitative  measurements  of  flow 
field  velocity  distribution  can  be  obtained.  In  cases  where 
solid  particles  are  used  as  tracer  materials,  errors  can 
result  when  tracer  particles  do  not  accurately  follow  the 
flow  streamlines. 

In  the  present  work,  the  use  of  a  relatively  new  flow 
visualization  and  measurement  technique  is  described.  The 
tracer  chosen  was  a  photochroraic  pyrospiran  dye  that,  while 
dissolved  in  a  hydrocarbon  solvent,  changed  in  color  from 
transparent  to  deep  violet  when  exposed  to  ultraviolet 
light.  A  pulsed  nitrogen  gas  laser  directed  at  a  fluid 
containing  the  dye  was  repeatedly  fired  into  the  fluid  pro¬ 
ducing  violet  traces  along  the  path  of  the  laser  beam. 
Subsequent  motion  of  the  fluid  deformed  the  traces  and  pro¬ 
duced  a  time  history  of  the  hydrodynamics. 


9 


Analysis  of  the  data  included  photographing  the  flow, 
digitizing  the  position  of  the  traces,  correcting  for  opti¬ 
cal  refraction,  and  computing  the  local  flow  velocities. 


A  high  powered  laser  was  constructed  based  on  a 
Blumenthal  circuit  [2]  producing  ultraviolet  light  at  a 
wavelength  of  3371  Angstroms.  The  internal  flow  of  dis¬ 
persed  droplets  of  a  hydrocarbon  solvent  growing  in  an 
immiscible  ambient  fluid  was  analyzed  to  verify  the  use  of 
the  technique. 
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II.  PREVIOUS  WORK 


Flow  visualization  through  the  use  of  photochromic  dye 
was  first  introduced  by  Popovich  and  Hummel  in  1966  [3]. 

They  found  that  by  dissolving  0.1%  by  weight  of  2(2,4- 
dinitro  benzyl) -pyridine  (DNBP)  in  951  ethanol,  blue  tracer 
lines  would  result  in  the  initially  clear  fluid  when  ex¬ 
posed  to  the  focused  beam  from  a  high  intensity  xenon  flash 
tube.  Quartz  plates  were  used  in  the  test  section  to  mini¬ 
mize  the  absorption  of  the  ultraviolet  light  from  the  flash 
tube.  The  pyridine  dye  shifted  in  color  from  transparent 
to  violet  within  3  microseconds  and,  if  not  reexposed,  re¬ 
verted  back  to  its  transparent  state  within  several  milli¬ 
seconds.  The  technique  was  used  to  look  at  the  viscous 
sublayer  in  turbulent  pipe  flow.  The  non- invasiveness  of 
the  technique  was  noted. 

The  production  of  DNBP  and  its  optical  characteristics 
are  explained  in  references  [4]  and  [5].  It  changes  color 
when  exposed  to  light  in  the  relatively  wide  range  of  ultra¬ 
violet  to  4000  Angstroms. 

Popovich  and  Hummel  [6]  employed  their  "flash  photoly¬ 
sis"  technique  using  pyridine  dye  in  ethanol  flowing 
through  a  smooth-walled  rectangular  aluminum  duct.  Photo¬ 
graphs  were  taken  of  1  mm  long  traces  produced  in  the 
ethanol  by  a  xenon  flash  tube  1.4  milliseconds  after  the 
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creation  of  the  trace.  The  deformation  of  the  traces  were 
analyzed  to  yield  the  velocities  in  the  laminar  sublayer 
and  in  the  buffer  layer.  It  was  noted  that  measurements 
beyond  Y+  =  17,  where 

Y+  =  Y  /T7p/v  (1) 

were  not  always  possible.  Results  were  compared  with 
Reichardt's  and  Diesler's  velocity  curves  and  showed  good 
agreement . 

Due  to  the  toxicity  of  DNBP,  the  use  of  1,3,3-tri- 
raethylindoline -6 1 -nitro-benzospiropyran  (TNBSP)  was  prefer¬ 
able.  This  dye  also  dissolved  easily  in  organic  solvents 
and  responded  to  ultraviolet  light  in  approximately  the 
same  manner  as  DNBP  after  ultraviolet  irradiation.  TNBSP 
transformed  from  a  transparent  state  to  a  violet  color 
(5250  Angstroms)  and  reverted  back  again  in  a  matter  of 
seconds.  Production  and  characteristics  of  TNBSP  are  ex¬ 
plained  by  Berman,  et  al  [7], 

Humphrey,  et  al  [8]  documented  the  fact  that  high 
fluxes  of  energy  from  an  ultraviolet  light  source  can 
drastically  alter  a  fluid  flow.  In  their  study,  a  giant 
pulse  ruby  laser  beam  fired  into  a  pendant  droplet  of 
chlorobenzene  containing  TNBSP  produced  violent  disturbances 
in  the  droplet.  Disturbances  were  attributed  to  tiny  gas 
bubbles  created  along  the  path  of  the  laser  beam  in  the 
fluid. 


Humphrey,  Hummel,  and  Smith  [11]  adapted  a  commercial 
multipulse  nitrogen  laser  to  replace  the  xenon  flashtube 
used  to  excite  the  photochromic  dye.  The  flow  visualiza¬ 
tion  technique  was  used  to  obtain  flow  patterns  inside  of 
growing  droplets  of  solvents  in  an  immiscible  ambient  fluid. 
Equations  were  given  to  facilitate  the  computation  of  veloc¬ 
ity  profiles  in  droplets  from  dye  traces.  Correcting  errors 
in  these  trace  positions  induced  by  optical  refraction  were 
discussed  in  reference  [9].  It  was  shown  that  errors  in 
trace  position  could  be  very  large  if  the  effects  of  refrac¬ 
tion  were  ignored.  Reference  [11]  contained  computed  veloc¬ 
ity  vector  distributions  along  the  central  core  of  a  slowly 
growing  droplet  of  an  organic  liquid  in  water.  The  sparse 
distribution  of  points  sampled  within  the  droplet  allowed 
approximate  pressure  profiles  to  be  computed. 

In  Culbreth,  Johnson,  and  Marschall  [12],  a  custom  built 
pulsed  nitrogen  laser  and  photochromic  dye  were  used  to 
study  droplets  of  mineral  spirits  rising  in  water  as  part 
of  a  study  of  direct-contact  heat  transfer.  High  speed 
movies  at  400  frames  per  second  were  made  of  developing 
droplets  as  the  pulsed  laser  produced  violet  traces  in  the 
droplet  interiors.  The  laser  design  was  used  in  the 
present  work  and  consisted  of  a  high-voltage  power  supply, 
a  source  of  pressurizing  nitrogen  gas,  a  vacuum  pump,  the 
laser  channel  made  up  of  plexiglas  box  with  copper 
electrodes  and  associated  electronics. 
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A  brief  history  of  the  pulsed  nitrogen  laser  used  for 
this  study  begins  with  the  development  of  an  inexpensive, 
powerful,  portable  device  by  James  Small  reported  in  refer¬ 
ence  [13].  His  simple  laser  design  involved  discharging 
two  copper  plates  across  a  1  cm  wide  gap  in  an  enclosure 
containing  nitrogen  gas  at  100  torr  of  pressure.  The  two 
plates  were  separated  by  several  sheets  of  mylar  plastic 
dielectric  above  a  ground  plate.  The  upper  plates  were 
charged  initially  to  approximately  20000  volts  and  one  plate 
was  periodically  grounded  producing  rapid  discharges  of 
electrons  through  the  enclosure  containing  the  nitrogen  gas. 
Collisons  between  electrons  excited  the  nitrogen  gas  into 
emitting  photons  with  a  wavelength  of  3371  Angstroms.  Pho¬ 
ton  emission  occurred  in  pulses  9  nanoseconds  in  duration 
during  which  50  to  100  kilowatts  of  ultraviolet  radiation 
were  released.  The  design  described  in  reference  [13]  re¬ 
lied  on  a  spark  gap  to  periodically  ground  out  one  of  the 
upper  plates  and,  subsequently,  fire  the  laser.  Commercial 
lasers  capable  of  producing  pulses  of  light  at  these  inten¬ 
sities  and  wavelength  cost  typically  on  the  order  of  $10000. 
In  reference  [2],  the  construction  of  a  small  pulsed  nitro¬ 
gen  gas  laser  is  described  costing  about  $30. 

Details  of  the  laser  constructed  for  the  present  study 
are  discussed  later.  Th.e  spark,  gap  used  to  fire  the  laser 
described  above  was  replaced  by  a  thyratron  tube  and  its 
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associated  electronics  as  described  by  Devlin  [10].  Devlin 
developed  a  liquid  cooled  pulsed  nitrogen  laser  for  non¬ 
linear  spectroscopic  measurements. 

The  history  of  the  pulsed  nitrogen-laser/photochromic 
dye  flow  visualization  technique  has  been  described.  The 
method  appeared  to  be  of  value  in  studying  a  wide  range  of 
fluid  dynamics  problems  where  non- intrusive  velocity  mea¬ 
surements  are  desired.  A  facility  has  been  developed  here 
at  NPS  as  a  part  of  this  work  using  this  measurement  tech¬ 
nique  and  the  strengths  and  weaknesses  of  the  technique  as 
a  measurement  tool  have  been  assessed.  Velocity  distribu¬ 
tions  have  been  obtained  in  growing  droplets  and  the  neces¬ 
sary  computer  codes  for  digitizing  and  reducing  data  are 
presented.  The  equations  necessary  for  correcting  optical 
distortion  and  for  computing  the  velocity  vectors  from  dye 
traces  are  given.  Construction,  operation,  and  safety  con¬ 
siderations  are  explained  for  the  pulsed  nitrogen  laser 
built  for  this  project. 


IS 


II 


EXPERIMENTAL  PROCEDURE 


The  apparatus  shown  in  Figure  1  consisted  principally  of 
a  vertical  3Mx3"xlO"  flat  wall,  plexiglas,  rectangular  par- 
allelpiped  containing  the  continuous  phase  (water).  The 
whole  of  the  experimental  equipment  was  constructed  of 
plexiglas  and  PVC  tubing.  Contamination  of  the  system  by 
dust  and  surface  active  agents  such  as  chemicals  and  grease 
was  carefully  avoided.  Droplets  of  mineral  spirits  were 
formed  from  a  single  polished  plexiglas  nozzle  immersed  in 
the  continuous  phase.  The  nozzle  was  pointed  so  that  the 
dispersed  phase  did  not  wet  the  nozzle  tip;  therefore,  the 
nozzle  diameter  could  be  considered  constant.  A  variable 
speed  positive  displacement  pump  was  used  to  circulate  the 
dispersed  phase  (mineral  spirits)  through  the  system.  The 
flow  rate  through  the  nozzle  was  regulated  by  a  combination 
of  pump  speed,  recirculation  valve  position,  and  cutout 
valve  position.  A  fluidic  capacitor  was  placed  in  line 
downstream  of  the  pump  and  upstream  of  the  flowmeter  as 
shown  in  Figure  2  to  dampen  the  pulses  that  were  generated 
by  the  pump  and  to  provide  a  storage  area  for  the  mineral 
spirits.  The  flow  rate  was  measured  using  a  previously 
calibrated  flowmeter.  The  mineral  spirits  was  recirculated 
through  the  system  until  equilibrium  temperature  was 
achieved  and  the  flow  rate  steadied  out  at  a  previously 
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selected  rate.  Upon  formation,  the  droplets  were  exposed 
to  the  ultraviolet  radiation  produced  by  the  nitrogen  laser. 
The  ultraviolet  light  from  the  laser  was  focused  through  a 
12  cm  focal  length  lens  until  only  a  very  thin  laser  line 
was  focused  in  the  droplet  just  above  the  nozzle  exit. 
Focusing  of  the  laser  beam  was  easily  accomplished  by  put¬ 
ting  a  few  drops  of  florishene  in  the  water  and  running  the 
laser  at  a  very  high  frequency.  This  allowed  a  visual 
check  on  both  the  focus  of  the  beam  and  the  exact  position 
that  the  beam  hit  the  droplet.  Upon  irradiation  of  the 
droplet,  the  trace  was  present  throughout  the  entire  history 
in  a  plane  containing  the  axis  of  symmetry  of  the  drop.  A 
high  resolution  television  camera  focused  through  a  10  power 
microscope  was  used  to  see  the  droplet  and  subsequent  laser 
traces  as  they  were  formed.  A  polaroid  picture  was  taken 
of  the  television  monitor  to  record  the  laser  traces.  Flow 
rates  were  varied  and  frequency  of  the  laser  was  varied  to 
obtain  the  best  overall  "picture''  of  the  forming  drop. 
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IV.  NITROGEN  LASER 


A.  DESIGN  CONSIDERATIONS 

This  thesis  required  the  use  of  a  high  powered  ultra- 
violet  nitrogen  pulsed  laser  to  cause  the  photochromic  dye 
to  change  color.  The  wavelength  of  the  laser  output  was 
337.1  am  with  peak  power  of  200- 300  megawatts  within  a  pulse 
width  of  approximately  9  nanoseconds.  The  pulsing  frequency 
was  variable  so  that  at  low  flow  rates,  the  droplets  would 
not  saturate  causing  the  laser  traces  in  the  droplet  to 
merge.  A  choice  of  a  high  voltage  thyratron  to  discharge 
a  capacitor  was  therefore  made. 

B.  LASER  PRINCIPLES 

The  basic  design  criteria  for  a  nitrogen  laser  is  to 
excite  as  many  molecules  as  possible  and  to  do  it  in  a 
short  period  of  time  (10  nanosec) .  The  closely  spaced  group 
of  laser  transitions  near  337.1  nanometers  in  the  nitrogen 
molecule  may  be  thought  of  as  arising  from  a  simple  four 
level  electronic  group.  See  Figure  3.  An  electron  in  the 
nitrogen  molecule  is  promoted  from  the  ground  state  to  the 
upper  laser  level  by  a  collision  with  a  free  electron  in  the 
high  voltage  discharge.  The  molecule  then  relaxes  to  the 
lower  laser  level  emitting  a  quantum  of  light  at  the  337.1 
nanometer  wavelength.  A  very  high  voltage  is  used  to  pump 
the  nitrogen  molecule  up  to  the  upper  laser  level.  The 
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laser  chosen  consisted  of  two  2S  cm  long  parallel  copper 
plates  separated  by  a  1  cm  gap.  A  rectangular  plexiglas 
enclosure  was  placed  around  the  plates  near  the  gap.  The 
enclosure  was  filled  with  nitrogen  gas  at  low  pressures  to 
provide  the  lasing  medium.  Glass  microscope  slides  with  a 
thickness  of  1.33  mm  were  attached  to  both  ends  of  the  en¬ 
closure  to  allow  the  laser  beam  to  pass  out  of  the  enclosure. 
The  laser  plates  were  separated  by  a  thin  mylar  sheet  with 
a  high  dielectric  constant  from  a  ground  plate.  The  plates 
were  configured  as  capacitors  by  an  electronic  circuit 
first  described  by  Blumenthal  [2],  By  pulsing  very  high 
voltages  across  the  gap  between  the  two  copper  plates  sur¬ 
rounded  by  nitrogen  gas,  the  gas  was  bombarded  by  electrons 
producing  the  plasma  which  emitted  a  quantum  of  light  in 
the  ultraviolet  range.  The  surface  area  of  the  copper 
plates  determined  the  magnitude  of  the  output  signal.  The 
larger  the  surface  area,  the  larger  the  capacitance  result¬ 
ing  in  a  larger  output  signal.  However,  increasing  the 
area  caused  the  time  to  charge  the  plates  to  increase. 

C.  LASER  CONSTRUCTION 
1 .  Laser  Channel 

The  copper  electrode  plates  of  the  capacitor  are  an 
integral  part  of  the  laser  channel.  Figure  4  shows  the 
dimensions  of  the  electrode  plates.  Plates  of  6.35  milli¬ 
meter  thick  plexiglas  were  used  to  construct  all  the  sides 
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of  the  laser  channel.  This  reduced  the  amount  of  ultraviolet 
light  emitted  to  the  surrounding  area  since  the  plexiglas 
absorbed  all  of  the  ultraviolet  light  produced  by  the  laser 
that  did  not  go  through  the  output  windows.  The  dimensions 
of  the  plexiglas  sides  are  shown  in  Figure  5.  Figure  6  shows 
how  the  plexiglas  pieces  and  the  electrode  plates  go  to¬ 
gether  to  form  the  laser  channel.  A  front  surface  mirror 
was  placed  at  the  back  of  the  laser  channel  to  double  the 
intensity  of  the  laser  output  by  reflecting  the  ultraviolet 
light  back  through  the  laser  channel  and  out  the  front 
window.  Silicon  sealer  was  used  to  seal  all  joints  and 
make  the  laser  channel  vacuum  tight. 

2.  Capacitor 

The  capacitor  was  made  from  thin  copper  plates  and 
layers  of  0.0075  inch  thick  mylar  plastic  sheets.  Figure  6 
shows  how  the  capacitor  was  assembled.  The  thyratron  lead 
was  soldered  to  the  top  capacitor  plate  and  a  ground  strap 
was  soldered  to  the  bottom  left  capacitor  plate.  The 
assembled  capacitor  was  then  bolted  together  using  phenolic 
particle  board  to  hold  the  copper  plates  and  mylar  sheets 
tightly  together.  The  measured  capacitance  was  found  to  be 
0.0045  microfarads.  Details  of  the  final  assembly  are 


shown  in  Figure  7. 
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3.  Electronics 

This  laser  design  utilized  an  EG$G  HY1102  grounded 
grid  thyratron.  Its  placement  in  the  laser  circuit  is  shown 
in  Figure  8.  All  electronic  components  were  installed  in  a 
cart  underneath  the  laser  channel  to  limit  the  lengths  of 
the  high  voltage  leads  to  the  laser  channel  and  to  make  the 
laser  semi-portable.  Figures  9  and  10  show  the  details  of 
the  triggering  circuits  for  the  thyratron. 

4.  Laser  Container 

The  laser  components  were  installed  in  the  cart  as 
shown  in  Figure  11.  This  cart  was  grounded  to  provide  the 
operator  protection  from  high  voltages  and  any  static 
electricity  buildup.  The  cart  also  attenuates  any  Radio 
Frequency  energy  generated  by  the  electronic  components. 

D.  AUXILIARY  SYSTEMS 

1.  Vacuum  System 

A  Welch  model  1402  DUO-SEAL  vacuum  pump  was  used  to 
maintain  the  laser  channel  between  the  ideal  26-28  inches  of 
vacuum.  The  suction  line  was  installed  at  the  rear  of  the 
laser  channel  and  a  0-3Q  inch  Hg  vacuum  gauge  was  placed 
between  the  laser  channel  and  the  vacuum  pump  close  to  the 
laser  channel. 

2.  Nitrogen  Gas  Supply 

Nitrogen  gas  was  supplied  to  the  front  of  the  laser 
channel  via  a  reducing  valve  and  cutout  valve.  Nitrogen 


pressure  to  the  cutout  valve  was  set  for  10  psi  and  the  cut¬ 
out  valve  was  used  to  regulate  the  laser  channel  pressure. 

3.  Power  Supplies 

The  thyratron  reservoir  filament  and  cathode  fila¬ 
ment  heaters  were  powered  by  a  115  v  AC  voltage  source 
through  a  transformer  which  reduced  the  voltage  to  6.2  volts. 
The  SCR  anode  voltage  was  supplied  by  an  HP  711A  0-300  VDC 
power  supply.  A  WAVETEC  pulse  generator  provided  the  +/- 
10  V  square  wave  signal  at  various  frequencies  to  the 
cathode  decoupler  network.  High  voltage  was  supplied  by  a 
HIPPOTRONICS  0-30  kV  power  supply.  Power  supply  connections 
are  shown  in  Figure  8. 

E.  LASER  OPERATION 

1.  Operational  Settings 

The  SCR  anode  voltage  was  set  to  250  volts.  The 
pulse  generator  was  at  the  desired  frequency  between  1  and 
30  Hertz.  Channel  pressure  varied  between  2-4  inches  abso¬ 
lute  of  mercury.  Alignment  of  the  output  mirror  was  done  by 
removing  the  mirror  and  sighting  down  the  laser  channel  to 
the  target.  The  mirror  was  installed  and  sighting  down  the 
laser  channel  in  the  opposite  direction  allowed  for  its 
alignment. 

2.  Firing  Sequence 

An  understanding  of  the  firing  sequency  is  essential 


to  the  safe  operation  of  the  electronic  circuitry. 


Initially  the  top  capacitor  plates  are  charged  to  a  voltage 
of  about  18,000  V  and  the  bottom  plate  is  grounded.  A  large 
negative  pulse  is  applied  to  the  thyratron  cathode.  This 
large  negative  pulse  causes  a  discharge  between  the  thyra¬ 
tron  grid  and  anode.  After  discharge,  the  top  left  capa¬ 
citor  plate  becomes  negatively  charged  with  respect  to  the 
top  right  capacitor  plate.  A  high  potential  between  the  two 
electrode  plates  will  result  in  breakdown  and  thus  excita¬ 
tion  of  the  nitrogen  molecule  by  electron  collisions.  An 
inductor  is  placed  in  line  between  the  left  and  right  capa¬ 
citor  plates  to  allow  the  charge  buildup  to  occur  on  the 
right  top  plate  and  to  allow  the  right  top  plate  to  ground 
completely  before  the  next  pulse  is  applied  to  the  left 
plate. 

3.  Turn-On  Procedure 

a)  Plug  in  the  power  cord  for  the  thyratron  heater 
elements.  A  10  minute  warmup  time  is  required  to  prevent 
damage  to  the  thyratron  when  operating. 

b)  Evacuate  the  laser  channel. 

c)  Connect  the  67.5  VDC  battery  to  the  terminals. 

d)  Adjust  the  nitrogen  gas  supply  to  maintain  26-28 
inches  Hg  vacuum  in  the  laser  channel. 

e)  Set  the  SCR  anode  voltage  for  250  Volts. 

f)  Turn  on  the  pulse  generator  and  set  to  the 
desired  frequency. 
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g)  Energize  the  high  voltage  power  supply.  Slowly 
raise  the  voltage  until  lasing  occurs  at  approximately  15  kV. 
Do  not  exceed  20  kV.  The  thyratron  can  be  damaged  if  the 
voltage  exceeds  20  kV.  NOTE:  The  high  voltage  will  in¬ 
crease  as  pulsing  frequency  decreases. 

4 .  Turn-Off  Procedure 

To  turn  off  the  laser  follow  the  turn-on  procedure 
in  reverse  order.  It  is  important  that  the  high  voltage  be 
lowered  before  lowering  the  frequency  because  you  can  ex¬ 
ceed  20  kV.  Turn  off  the  nitrogen  supply  and  let  the  pres¬ 
sure  bleed  down  before  turning  off  the  vacuum  pump  to 
eliminate  the  possibility  of  creating  a  pressure  inside  the 
laser  channel  greater  than  atmospheric.  This  situation  can 
lead  to  pinhole  leaks  forming  in  the  laser  channel. 

F.  LASER  SAFETY 

Ultraviolet  radiation  can  cause  serious  damage  to  the 
operator's  eyes.  To  reduce  the  hazard,  the  laser  was  set 
up  in  a  controlled  access  area.  Only  those  people  wearing 
plastic  safety  glasses  were  allowed  in  the  room  while  the 
laser  was  in  operation. 

Since  very  high  voltages  were  used,  the  laser  container 
was  grounded  to  prevent  a  shock  hazard.  It  also  prevented 
a  static  electricity  buildup  on  any  components. 
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V.  DATA  ANALYSIS 


The  photograph  showing  the  clearest  laser  traces  was 
digitized  using  an  HP-9845  computer  system.  Figure  12  shows 
the  data  points  used  in  the  analysis. 

The  first  step  was  to  correct  all  the  data  points  for 
the  difference  in  refraction  between  the  mineral  spirits 
and  water.  The  correction  is  necessary  because  a  curved 
surface  separating  two  medias  of  different  indices  of  re¬ 
fraction  will  cause  the  apparent  displacement  of  any  point 
within  one  or  the  other  mediums  when  viewed  through  the 
interface. 

A.  ASSUMPTIONS 

The  following  assumptions  are  made: 

1)  The  drop  is  symmetrical  about  its  vertical  axis  and  its 
profile  can  be  accurately  described  by  the  empirical  equa¬ 
tion  given  by  Poutanen  and  Johnson  [2]  for  growing  bubbles. 

Rn9  =  1  (2) 

where:  R,9  =  polar  coordinates  of  the  profile 

n  *  arbitrary  parameter  chosen  to  fit  any 
particular  droplet 

Figure  13  shows  the  drop  geometry  for  Equation  (2). 


2)  The  point  object  is  located  in  a  plane  containing  the 
drop  axis;  the  said  plane  being  perpendicular  to  the  normal 
of  the  camera  lens. 

3)  All  rays  of  light  proceeding  from  point  objects  within 
the  drop,  upon  crossing  the  interface  can  be  assumed  to 
follow  paths  parallel  to  the  normal  of  the  camera  lens. 

4)  The  center  of  gravity  of  the  growing  drop  and  the  camera 
lens  central  axis  remain  at  the  same  height  during  drop 
growth. 


B.  CALCULATION  OF  OPTICAL  DISPLACEMENTS 

The  direct  application  of  Equation  (2}  in  polar  form  is 
very  inconvenient  for  the  present  analysis.  A  transforma¬ 
tion  into  the  Z -Y  coordinate  system  was  made  taking  the 
apex  of  the  drop  as  the  origin.  The  result  of  the  trans¬ 
formation  is: 

(Y2  ♦  Z  2)ARCTAN(^-)  =  ln  (3) 

l  lx 

where:  Z.  =  Z  -  R  (4) 

1  o 

1  -  Dm/d’  (5) 

R  *  distance  from  the  polar  coordinate  origin 
to  the  new  origin 

n,d'  =  empirical  constants  determined  as  in  [9] 
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It  can  be  shown  that  the  optical  corrections  can  be  expressed 


as  : 

AZ  =  TAN(r-i)  /C7 -R2  VTAN (r)  (6) 

AY  *  (TAN(r-i) /TAN(r) ) x  R2  (7) 

where:  r  -  ARCOS(C/(R1  STTUTm^ 77  (8) 

i  =  arcsin(nr/n£  SIN(r))  (9) 

C  *  /R~2  -  R., 2  (10) 

R^  ■  Distance  from  center  to  edge  of  bubble. 


R2  =  Distance  from  center  to  edge  of  point  to  be  corrected. 

n  *  Index  of  refraction  of  water 
r 

*  Index  of  refraction  of  mineral  spirits 

N+2 

(Zf-R0)  (R^+CZf-R^2)^-  cos2(JlN(R12+(Zf-R0)2)  ’*/2)+rn\  (Zf-RQ) 
“l  =  S+2 

R1(R12+  CZ£-R0)2)"7"cos2CiN(R12+  (Z£-R0)2)'N/2)  -  nilNCZf-Ro)3 

(11) 

Z^  =  Vertical  coordinate  of  point  F 
The  optically  corrected  data  points  are  then  determined  by: 


X  +  AY 

(12) 

Z  +  AZ 

(13) 

Figure  14  shows  the  position  of  the  corrected  data  points. 
Figure  IS  shows  the  position  of  the  corrected  and  uncorrected 
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data  points.  Figure  18  is  the  HP-9845  computer  program 
that  was  used  to  solve  the  above  equations  for  the  optical 
corrections . 

C.  VELOCITY  VECTOR  DETERMINATIONS 

Due  to  symmetry,  only  the  right  half  of  a  droplet  need 
be  considered  as  shown  in  Figure  16.  The  Y  axis  lies  along 
the  center  line  of  the  nozzle  and  is  given  the  value  0  at 
the  nozzle  exit,  Y  is  positive  in  the  direction  of  flow. 

The  radial  coordinate  is  labeled  X,  and  is  positive  away 
from  the  centerline. 

Quantitative  measurements  of  fluid  flow  inside  the  drop¬ 
let  were  limited  to  a  core  along  the  Y  axis  with  the  same 
diameter  as  the  inside  diameter  of  the  nozzle.  In  this  well 
defined  region,  the  shape  of  any  induced  trace  i,  after  cor¬ 
rection  for  optical  distortion,  was  found  to  be  excellently 
represented  by  a  polynomial  of  degree  3  i.e., 

F (X)  *  a(I)X3  +  b(I)X2  +  c(I)X  +  d(I)  Cl*) 

where  a(I),  b(I),  c(I),  and  d(I)  are  time  dependent  and  are 
calculated  hy  a  least  squares  fit  of  Equation  (14)  to  the 
experimental  values. 

The  fact  that  flow  within  the  growing  droplet  is  not 
undirect ional  except  at  very  high  flow  rates  complicates 
the  calculation  of  fluid  element  velocities  from  the  induced 
traces.  Two  traces  within  the  drop  at  time  t  are  shown  as 
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dotted  lines  in  Figure  16.  At  a  later  time  t+At,  the  same 
traces  are  shown  as  continuous  lines  .  It  is  possible  to 
follow  the  change  of  position  of  an  element  of  fluid  from 
point  A  to  point  B  using  a  mass  balance  technique  as  follows. 
Consider  the  displacement  of  an  element  of  fluid  (PI,  P2, 

P3,  P4)  contained  between  traces  1  and  2  during  an  interval 
of  time  At  as  shown  in  Figure  16.  As  the  element  moves  for¬ 
ward  to  its  new  position  it  is  distorted,  some  of  the  fluid 
from  the  moving  element  being  squeezed  out  in  the  radial 
direction.  Provided  the  fluid  element  is  small,  an  approxi¬ 
mation  to  its  vector  velocity  is: 

V  =  XI/At  (15) 

and  the  direction  of  its  movement  can  be  characterized  by 
the  angle  oc. 

Given  trace  1  and  2  , 

Let  A^  =  X  coordinate  of  point  A 
Let  Ay  *  Y  coordinate  of  point  A 

Then 

Ay  *  l/2[(a(2)-a(l))Ax3+(b(2)-b(l))Ax2 

+  (c(2)-c(l))Ax  +  d(2)-d(l)]  +  F(Ax)  (16) 

Though  distorted  due  to  a  combined  radial  and  axial  motion 
of  the  fluid,  the  volume  of  the  fluid  element  remains  the 
same  during  an  interval  of  time  At.  Since  the  trace  con¬ 
taining  the  distorted  element  can  also  be  fitted  to 
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equations  of  the  type  given  by  Equation  (14) ,  where  the 
constants  are  now  found  for  that  element  of  volume,  the 
volumes  are  set  equal  to  each  other  and  the  polynomial  is 
solved  for  the  X  coordinate  of  point  B.  Therefore  from  the 
theorem  of  PAPPUS: 

Let  ■  X  coordinate  of  point  B 
Let  By  =»  Y  coordinate  of  point  B 

V(l,2)  3  Ca(2)-aCl))Ax5*Cb(2)-b(l))Ax4 

+  (c(2)-c(l))Ax3+(d(2)-d(l))Ax2  (17) 

V  (3 ,4)  =  CaC4)-aC3))3x5  +  CbC4)-b(3))Bx4 

+  (c(4)-c(3))Bx3  +  Cd(4)-dC3))Bx2  (18) 

Now,  since  V(l,2)  3  V(3,4)  Equation  CIS)  can  be  solved  for 
Bx-  Once  this  value  is  found,  substitution  into  the  follow¬ 
ing  equation  gives  By 

3y  3  l/2[(a(4)-a(3))Bx3  +  (b  (4) -b  (  3)  )  B^ 

+  CcC4)-cC3))Bx+  (d(4)-d(3))]+F(Bx)  (19) 

The  velocity  can  now  be  expressed  as: 

V  3  /(Bx-Ax)2  +  (By-Ay)2  (20) 

and  the  angle  a  can  be  found  by: 

a  3  ARCTAN((By-Ay)/(Bx-Ax))  (21) 
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Figure  19  is  the  HP-984S  computer  program  that  was  used  to 
solve  the  above  equations  for  the  velocity  vectors  and 
angles . 


31 


VI .  RESULTS 


Figure  17  shows  the  velocity  profiles  for  the  4  traces 
that  were  obtained  during  this  procedure.  As  can  be  seen, 
the  velocity  decreased  as  the  fluid  moved  away  from  the 
nozzle  in  the  vertical  direction;  and  the  velocity  along  a 
given  trace  also  decreases  as  you  move  away  from  the  center 
line.  The  magnitude  of  the  peak  velocity  coming  out  the 
nozzle  was  computed  to  be  97.6  micrometers  per  second.  This 
corresponds  to  a  Reynolds  number  of  approximately  0.05. 

The  velocity  profile  looked  very  much  like  the  laminar 
velocity  profile  at  the  nozzle  exit  which  is  parabolic  in 
shape.  Table  I  shows  the  coordinates  and  magnitude  and 
direction  of  the  velocity  vectors  using  laser  traces  1-2, 
2-3.  The  origin  of  the  coordinate  system  is  the  center  of 
the  nozzle  at  the  nozzle  exit.  The  magnitude  of  the  veloc¬ 
ity  is  unitless  and  is  equal  to  the  computed  velocity 
divided  by  the  peak  velocity  in  the  nozzle.  Likewise, 

Table  II  shows  the  results  using  laser  traces  2-3,  3-4  and 
Tahle  III  shows  the  results  using  laser  traces  1-3,  2-4. 
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A.  ERROR  ANALYSIS 

All  positional  measurements  in  the  droplets  were  based 
on  the  observed  nozzle  bore  diameter.  The  accuracy  of  the 
bore  diameter  was  approximately  2.7%.  The  droplet  profiles 
and  position  of  points  on  traces  were  digitized  from  a 
Polaroid  picture  of  droplets  obtained  through  a  videoscreen. 
The  effects  of  parallax  were  minimized  by  placing  the  camera 
as  far  from  the  inlet  nozzle  as  possible.  Humphrey  and 
Hummel  [9]  reported  that  parallax  error  is  negligible  if 
the  droplet  traces  are  observed  from  more  than  25  cm  and 
readings  are  confined  to  within  3  mm  of  the  center  of 
gravity  of  the  droplet.  In  this  study,  the  camera  was 
placed  approximately  15  cm  away  from  the  droplet  and  traces 
were  viewed  within  1  ram  of  the  droplet  center.  The  error 
caused  by  the  nozzle  bore  diameter  far  exceeded  the  digitiz¬ 
ing  error  of  0.008  mm. 

The  error  induced  in  determining  each  point  on  a  trace 
in  (X,Y)  coordinates  was  influenced  by  errors  in  maximum 
droplet  diameter,  point  coordinates  R2  and  Zp  described 
earlier,  and  the  indices  of  refraction  of  mineral  spirits 
and  water. 

The  errors  in  indices  of  refraction  were  assumed  to  be 
small  and  in  the  correction  for  optical  refraction,  posi¬ 
tional  errors  predominated.  An  analysis  was  done  to  compute 
the  error  in  position  for  points  on  the  traces  in  the 


33 


droplets  after  correcting  for  refraction.  The  uncertainty 
analysis  technique  detailed  by  Holman  [14]  was  used.  The 
points  obtained  from  photochromic  dye  traces  that  were  ed 
to  compute  velocities  were  all  confined  within  an  area  no 
wider  than  the  inlet  nozzle  diameter  and  far  from  the  a.  ex 
of  the  droplet.  Within  this  area,  point  positions  were 
known  within  +_  0.005  mm  after  correcting  for  refraction. 

The  corrected  traces  were  curve  fitted  using  a  third 
order  polynomial.  The  fit  was  very  good  with  a  standard 
deviation  of  under  0.001  mm.  The  time  interval  between 
firings  was  accurately  controlled  by  a  crystal  based  func¬ 
tion  generator  with  an  accuracy  of  about  +  0.0001  seconds. 
All  computations  on  the  HP-9845  were  done  in  double  preci¬ 
sion  arithmetic  to  minimize  errors. 

The  estimated  accuracy  of  the  fluid  velocities  was,  on 
the  average,  0.07  mm/s  or  6.3%  of  the  average  fluid  velocity 
at  the  nozzle  tip. 

It  should  be  noted  that  errors  will  rise  very  rapidly 
if  data  is  taken  near  the  periphery  of  the  droplet  due  to 
the  traces  appearing  to  be  very  close  to  one  another. 
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VII.  CONCLUSIONS 


Velocities  were  computed  using  a  nitrogen  laser/flow 
visualization  technique.  The  technique  was  used  to  investi¬ 
gate  the  velocity  profile  in  the  core  of  a  slowly  forming 
droplet  (Re  *  Q.Q5).  Velocities  along  the  centerline  of 
the  droplet  were  seen  to  decrease  from  a  peak  velocity  in 
the  nozzle  of  2.22  ram/s  to  95.61,  87.8%,  and  68.11  of  the 
peak  velocity  as  one  moves  toward  the  droplet  apex.  Peak 
velocity  in  the  nozzle  equaled  twice  the  average  velocity 
as  determined  from  parabolic  laminar  profile  in  a  circular 
pipe.  Computed  velocities  agreed  well  with  expected  veloc¬ 
ities  and  with,  the  restrictions  of  laminar  flow  and  small 
flow  geometry,  this  measurement  technique  appears  to  be 
accurate  and  useful. 


VIII.  RECOMMENDATIONS 

Since  a  "home-made"  laser  was  used  in  this  experiment, 
only  simple  hydrodynamic  problems  can  be  solved.  The  ab¬ 
sorption  of  the  laser  light  by  the  mineral  spirits  necessi¬ 
tated  the  use  of  a  very  small  diameter  nozzle.  Only  with 
high  powered  commercial  grade  lasers  can  this  technique  be 
expanded  to  more  challenging  hydrodynrmic  problems.  Since 
only  polaroid  pictures  were  taken  the  rate  of  droplet  growth 
was  severely  limited.  This,  in  turn,  limited  the  laser 
frequency  so  that  saturation  did  not  occur. 
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Figure  6 
Laser  Channel 
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Figure  8 

Thyratron  Placement  in  Laser  Circuit 


Decoupler  Network 


Triggering  Circuit 


Figure  11 


Electronic  Component  Configuration 


Figure  13 


Geometry  for  Equation 
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Figure  14 

Corrected  Data  Points 
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Figure  15 

Optically  Corrected  and  Uncorrected  Data  Points 
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Figure  16 


Velocity  Vector  Computations 
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Figure  17 

Velocity  Profile  inside  the  Droplet 
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10  OPTION  BNSE  1 
1  1  DEG 

20  DIM  R2C  100  <  ,  PI  ‘  100)  ,  2»  '  lOO  > ,  1O0.>  ,  100  • 

30  NS SIGN  #1  TO  “ BUB  1 " 

40  PEND  til ,  1 
50  FOR  1*1  TO  24 
60  REND  ii;X<I},Y(I> 

70  R2U>*X<I> 

80  Zf  <I>*4039-YU> 

90  NEXT  I 

91  NSSIGN  #1  TO  “  R 1 1  :F8,0*' 

92  REND  il,  1 

100  FOR  1-1  TO  24 
110  REND  #l;Rl<I'- 
120  NEXT  I 
130  L*20 1 0 

140  N*5 . 0 

150  R0* 1 608 . 9 

160  £tai*1.48 

170  Et  ar* 1.41 
180  D«1ta=L'N 

190  CRENTE  "BUBLEl : F3" , 1 , 16*50 
200  NSSIGN  ttl  TO  “BUBLEl “ 

210  REND  til,  1 

220  FOP  1-1  TO  24 

230  C*SQR<R1<  I  > -“>2 -R2 <  I  }  s  2  > 

240  N 1 pha*ZF< I } -R0 

250  B«t  a*R  1  <  I )  '2+N 1  pha-'Z 

251  Q*Nlph**B*ta-'3. 5 

260  M0-Q*COS < D«  1  2>  >  '2 

270  M0*M0+N*De 1 1  a*R 1 < I >  *N 1 pha  '2 

271  Q1»R1<I)»B*U"3.3 

280  M2=Q  1  *COS  <  D«  1 1  a*Be  t  a  '  <.  -N  -  2  >  )  ,v2 
290  M2-M2-N*D«lt**N1pha'3 

300  M 1 -M0' M2 

310  Rh*t*NCS<.C-'<Rl  <  I  >*SQR<  l*MlM-2>  >  >  ) 

320  Ih*t=NSN<Et*r'Et ai *SIN<Rhat  >  > 

330  D«ltay*TNN<Rh4t-Ihat }^TNN<Rhit  >*R2< I > 

340  D*  1  t  az*TNN  <  Rhat - 1  hat  >*SQRCC"2-R2<  I  )  A2''TNN  ( Rhat 
350  X<I)*X(I)-D«lt»y 
360  YU)»Y(I)-D«lta2 
370  PRINT  #1 ; X< I ) , Y<  I  ) 

380  PRINT  X<I),D*Uay,Y<I>,D*Uaz 

390  NEXT  I 

400  DISP  “COMPLETE" 

410  STOP 
428  END 


Figure  18 

Program  to  Determine  Optical  Corrections 
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le  fll— 9. 0056E-7 
20  R2— 1 . 00831E-8 

30  R3— 5.3438478E-7 

40  R4*9. 857368E-8 

30  Bl». 00005203 
<0  82*. 80031211875 

70  83— .00027479 

80  84*-. 001 15880083 

90  Cl  — .  17388898 
100  C2— .23133398 

110  C3— .05371831 

120  C4-. 23938038 

130  01*1971 . 8236 

140  02*2423.3385 

130  03*2718. 73929 

180  04*2943.295034 

170  21*02-01 

180  22*82-81 
190  Z3*C2-C1 

200  24-02-01 

210  23*03-02 

220  26*83-82 

230  27*C3-C2 

240  28*03-02 

230  CRE0TE  "BBVEL1 “ , 1 , 2830 
260  ASSIGN  «1  TO  “BBVEL1  ” 

270  RERO  *1.1 
280  L2*0 

290  L1*0 

291  DEG 

292  BUFFER  *1 
300  FOR  1-1  TO  89 

310  Ry*<21»LlA3*Z2*Ll  A2  +  Z3*Ll*Z4>^2+fll*LlA3  +  Bl#LlA2+Cl*Ll-t-Dl 

320  Co*f*Zl*LlA5*Z2*HA4*Z3*LlA3+Z4*LlA2 
330  FOR  J-l  TO  2000 

340  Ck»Z3*LZA5'*’ZS*L2''4+Z7*L2A3*28*L2A2-Co*f 

330  IF  Ck  >0  THEN  GOTO  380 
360  L2-L2*l 

370  NEXT  J 

300  By*<Z3*L2A3*Z6*L2A2*Z7*L2+Z8>.'2'*A2*L2A3  +  B2*L2A2+C2*C2't-D2 

390  Ab*. 23*SQR<  <L2-L 1 )  A2*<  Bv-Ay > A2  > 

400  IF  L2-L1-0  THEN  GOTO  430 
410  Th*t **ATN<  <  By-Ay ) / <  L2-L1 ) ) 

420  GOTO  440 
430  Th*t**98 

440  PRINT  Ll.Ay.L2.By,  Rb,Th*t* 

430  PRINT  *i;Ll,Ay,L2,8y 

451  IF  L2<200  THEN  L2*0 

432  IF  L2>200  THEN  L2-L1-20O 

460  Ll-Ll+10 

470  NEXT  I 

480  BISP  “COMPLETE" 

490  RSSIGN  *1  TO  • 

300  STOP 
310  ENO 


Figure  19 

Program  to  Determine  Velocity  Vectors 
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TRACES  1 -2-3 


X 

Y 

VELOCITY 

THETfl 

30.000 

2198.592 

.956 

89.846 

60.000 

2192.351 

.957 

88.923 

90.000 

2186.294 

.958 

87.852 

120.000 

2180.235 

.957 

87.087 

1S0.000 

2174.050 

.956 

86. 165 

180.000 

2167.574 

.953 

85.235 

210.000 

2160.653 

.949 

84.450 

240.000 

2153.131 

.943 

83.492 

270.000 

2144.855 

.935 

82.330 

300.000 

2135.669 

.924 

81.172 

330.000 

2125.420 

.910 

79.778 

360.000 

2113.952 

.894 

78.312 

390.000 

2101.110 

.873 

76.573 

420.000 

2086.741 

.849 

74.508 

490.000 

2070.689 

.821 

72.247 

480.000 

2052.800 

.790 

89.733 

310.000 

2032.920 

.756 

66.893 

540.000 

2010.892 

.720 

63.631 

570.000 

1986.564 

.684 

60. 124 

600.000 

1959.780 

.648 

56.327 

630.000 

*  1930.386 

.613 

32 . 204 

660.000 

1898.227 

.580 

47.729 

690.000 

1863.148 

.532 

42.903 

720.000 

1824.993 

.527 

38.309 

750.000 

1783.613 

.504 

34.214 

780.000 

1738.848 

.484 

30.261 

810.000 

1690.344 

.466 

26.620 

840.000 

1638.548 

.450 

24.287 

870.000 

1582.704 

.435 

22.807 

- 

Table  I.  Numerical  results  of  velocity  vector  computations 
for  laser  traces  1-2,  2-3. 
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TRRCES  2-3-4 


30.000 
60.000 
90.000 
120.000 
ISO. 000 
180.000 
210.000 
240.000 
270.000 
300.000 
330.000 
360.000 
390.000 
420.000 
450.000 
480.000 
S10.000 
540.000 
570.000 
800.000 
630 
660 
690 
720.000 
750.000 
780.000 
810.000 
840.000 
870.00** 


2571.074 
2566.494 
2561.820 
2556.927 
2551.508 
2545.876 
2539.655 
2532.527 
2524.738 
2515.868 
2505.894 
2494.687 
2482.120 
2468.069 
2452.405 
2435.003 
2415.735 
2394.475 
2371.09 7 
2345.473 
2317.478 
2286.985 
2253.867 
2217.997 
2179.249 
2137.497 
2092.613 
2044.471 
1992.945 


VELOCITY 

THETR 

.681 

90.000 

.695 

89.138 

.705 

88.521 

.712 

87.917 

.715 

87.318 

.715  i 

86.714 

.713 

86.099 

.707 

85.467 

.699 

84.598 

.688 

83.690 

.674 

82.508 

.657 

81.250 

.637 

79.899 

.615 

77.921 

.589 

76.020 

.563 

73.322 

.536 

70.601 

.508 

67.515 

.482 

63.989 

.458 

60.376 

.437 

56.734 

.418 

53.686 

.403 

50.919 

.390 

48.587 

.381 

46.864 

.375 

46.646 

47.452 

48.630 

.383 

50.953 

Table  II.  Numerical  results  of  velocity  vector  computations 
for  laser  traces  2-3,  3-4. 


TRACES 

1 -3-2-4 

X 

Y 

VELOCITY 

THETR 

30.009 

2344.307 

.876 

89.831 

60.000 

2340.773 

.883 

88.997 

90.000 

2336.921 

.  887 

88.172 

120.000 

2332. S33 

.888 

87.517 

150.000 

2327.790 

.888 

86.693 

.  1B0.000 

2322.276 

.885 

86.030 

210.000 

2315.971 

.881 

85.187 

240.000 

2308.760 

.875 

84.327 

270.000 

2300.523 

.866 

83.442 

300.000 

2291.142 

.855 

82.349 

330.000 

2280.500 

.841 

81.209 

360.000 

2268.490 

.823 

79.821 

390.000 

2254.962 

.803 

78.150 

420.000 

2239.830 

.780 

76.353 

450.000 

2222.966 

.754 

74. 179 

480.000 

2204.250 

.724 

71.779 

510.000 

2183.567 

.693 

68.838 

540.000 

2160.798 

.662 

65.778 

570.000 

2135.824 

.632 

62.588 

600.000 

2108.529 

.602 

59.266 

630.000 

2078.794 

.574 

55.825 

660.000 

2046.502 

.549 

52.294 

690.000 

2011.535 

.526 

49.203 

720.000 

1973.774 

.507 

46.250 

750.000 

1933. 102 

.490 

44. 1 16 

780.000 

1899.401 

.475 

42.459 

810.000 

1842.554 

.463 

41.429 

840.000 

1792.441 

.456 

41.182 

870.000 

1738.947 

.451 

42.590 

1 

_ _ _ 1 

Table  III.  Numerical  results  of  velocity  vector  computati 
for  laser  traces  1-3,'  2-4. 
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